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Based on the Skyrme energy density functional, the spatial distribution of the symmetry energy of a finite nu- 
cleus is derived in order to examine whether the symmetry energy of a finite nucleus originates from its interior 
or from its surface. It is found that the surface part of a heavy nucleus contributes dominantly to its symmetry 
energy compared to its inner part. The symmetry energy coefficient a sym (A) is then directly extracted and the 
ratio of the surface symmetry coefficient to the volume symmetry coefficient k is estimated. Meanwhile, with 
the help of experimental alpha decay energies, a macroscopic method is developed to determine the symme- 
try energy coefficient of heavy nuclei. The resultant a sym (A) is used to analyze the density dependence of the 
symmetry energy coefficient of nuclear matter around the saturation density, and furthermore, the neutron skin 
thickness of 208 Pb is deduced which is consistent with the pygmy dipole resonance analysis. In addition, it is 
shown that the ratio k obtained from the macroscopic method is in agreement with that from the Skyrme en- 
ergy density functional. Thus the two completely different approaches may validate each other to achieve more 
compelling results. 

PACS numbers: 21.65.Ef, 21.65.Cd, 21.60Jz, 23.60.+e 



I. INTRODUCTION 

Great attention has been paid to the nuclear equation of state 
(EOS) of isospin asymmetric nuclear matter, in particular 
the nuclear matter symmetry energy coefficient S(p). The 
density-dependent symmetry energy coefficient plays a cru- 
cial role in understanding a variety of issues in nuclear physics 
as well as astrophysics, such as the heavy ion reactions [1-5], 
the stability of superheavy nuclei [6 , the structures, compo- 
sition and cooling of neutron stars [|7 -U0| an d even some new 
physics beyond the standard model 1211 . The energy per 
particle in nuclear matter is e(p, 6) = e(p, 0) + S (p)5 2 + 0(5 4 ) 
with the density p — p„ + p p and isospin asymmetry 6 = 
(fin - Pp)lp- The density-dependent symmetry energy co- 
efficient S(p) can be expanded to second order by S(p) = 
So + Le + K sym e 2 /2 with e = (p - po)/3po, where So is the 
symmetry energy coefficient at the nuclear matter saturation 
density, L = 3pdS (p)/dp\ Po and K sym = 9p 2 d 2 S /dp 2 \ Po are the 
slope and curvature parameters governing the density depen- 
dence of S(p) around the saturation density p<). 

Various independent investigations have been carried out to 
constrain the density dependence of the nuclear matter sym- 
metry energy coefficient. A microscopic Brueckner-Hartree- 
Fock approach using the realistic Argonne V18 nucleon- 
nucleon potential plus a phenomenological three-body force 
gives a value of the slope parameter L = 66.5 MeV lfl3ll . 
In a description of isospin diffusion data using a Boltzmann- 
Uehling-Uhlenbeck (BUU) transport model (IBUU04) with 
free-space experimental nucleon-nucleon (NN) cross section, 
Chen et al. extracted a nuclear symmetry energy coefficient 
of S(p) ~ SoiplpoY with y — 1.05 by comparing the the- 
oretical results with the experimental data II 1411 . Including 
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the medium-dependent NN cross section, the isospin diffu- 
sion data lead to a softer symmetry energy coefficient with 
y = 0.69. The pygmy dipole resonance (PDR) of 208 Pb 
analyzed with Skyrme interactions gives L = 64.8 + 15.7 
MeV lfl5ll . Liu et al. extracted the symmetry energy coeffi- 
cients for the finite nuclei with mass number A — 20 - 250 
from more than 2000 measured nuclear masses [16]. With 
the semiempirical relationship between the symmetry energy 
coefficients a sym (A) of finite nuclei and S(p) of nuclear mat- 
ter, they obtained the slope parameter L - 53 - 79 MeV 
at the normal density. The analysis of isospin diffusion and 
double ratio data involving neutron and proton spectra by 
an improved quantum molecular dynamics transport model 
suggests S(p) = 12.5 (p/p () ) 2/3 + C p (p/p Y filwith 7 = 
0.4-1 .05. Danielewicz and Lee obtained the L values ranging 
from 78 to 111 MeV [18]. Centelles et al. recently reported 
the values of L = 75 + 25 MeV and later revised to a narrow 
window of L — 45 - 75 from neutron skin thickness of nuclei 
lfl9l l20ll . With a correlation for the symmetry energy at the 
saturation density So, the slope parameter L and the curvature 
parameter K sym based on widely different mean-field interac- 
tions, Ref. [21] suggested that the L value is 56 + 24 MeV. 
The neutron skin thickness AR„ p given by the difference of 
neutron and proton root-mea n-sq uare radii, correlates linearly 
with the slope parameter L Il22l424ll . Therefore, a measure- 
ment of AR„ p with a high accuracy should be a strong con- 
straint on the density dependence of S(p). Nevertheless, the 
density dependence of S (p) at subnormal density is currently 
still an open question, and further investigations are required. 

Recently, the symmetry energy of finite nuclei has been 
widely investigated because with the help of it one may gain 
some information on the density dependence of S(p). The 
most commonly used methods are the liquid drop models 
combined with nuclear masses and the leptodermous expan- 
sion based on the self-consistent mean-field theory Ilia, |25| - 
l27ll . In this study, we extract directly the symmetry energy 
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coefficient of finite nuclei a sym (A) by using the Skyrme energy 
density functional and mainly focus on the spatial distribution 
of the symmetry energy of a heavy nucleus. We would point 
out that compared to previous studies, we employ a different 
approach that the leptodermous expansion is not applied to a 
mean-field theory. Another purpose of the present work is to 
extract a sym (A) based on experimental alpha decay energies of 
heavy nuclei and investigate the density dependence of S(p). 

This work is organized as follows. In Sec. II, the approach 
to construct the symmetry energy density functional of finite 
nuclei is presented. a sym (A) and its spatial distribution is cal- 
culated and analyzed. In Sec. Ill, a sym (A) is extracted by using 
the experimental alpha-decay energies and the density depen- 
dence of S (p) is then determined. Finally a summary is given 
in Sec. IV. 



The fo-?3, X0-X3 and a are Skyrme interaction parameters re- 
sulting from a fit to the binding energies and charge radii of a 

number of known nuclei. p q , T q and J „ are the local nucleon 

density, kinetic energy density and spin density for neutron or 

— > 

proton respectively, whereas p, r and J are the corresponding 
total densities. After a self-consistent SHF calculation, one 
obtains the proton and neutron distributions p p (r) and p„(r) 
in a nucleus. By using the obtained SHF density distribu- 
tions, the total energy of the nucleus can be calculated. In 
order to avoid the implicit use of single particle wave func- 
tions for calculating the contributions of the kinetic energy to 
the symmetry energy, we express the kinetic energy density as 
a function of the density p q {r) by using the extended-Thomas- 
Fermi approximation to second order terms for reasons of be- 
ing enough for numerical convergence (HH, and it is written 
as (excluding the spin part) 133143511 



II. SYMMETRY ENERGY COEFFICIENT OF FINITE 
NUCLEI FROM THE SKYRME ENERGY DENSITY 
FUNCTIONAL 

The Skryme interaction has been a powerful tool in the in- 
vestigation of finite nuclei and nuclear matter [29-32]. In the 
standard Skyrme-Hartree-Fock (SHF) method, the interaction 
is a zero-range, density- and momentum-dependent form. The 
expression for the total energy density functional is 
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where q denotes neutron or proton. Here, we stress that the 
extend-Thomas-Fermi approximation is only adopted for cal- 
culating the total kinetic energy after the self-consistent SHF 
solution. Let p„ = p p = p/2 in the above formula (symmetry 
case), the density functional is given by < H^- Therefore, ex- 
cluding the Coulomb energy and the spin energy in 'H - "Ho, 
the density functional for the symmetry energy is written as 
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In a self-consistent SHF calculation, one obtains the nucleon 
density distribution finally. With the nucleon density distri- 
butions, the Ht, Hv an d 'Hgrad can be estimated accordingly. 
The total symmetry energy of a finite nucleus is defined as 
£ sym (A) = a sym (A)(N - Zf/A = f Q °°?{ sym dV. Here the sym- 
metry energy coefficient a sym (A) includes the volume and sur- 
face symmetry energy coefficient. 

Indeed, in order to extract the nuclear symmetry energy that 
relates solely to the nuclear force, one should subtract the ef- 
fect due to the fact that the Coulomb interaction effectively 
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FIG. 1: (Color Online) Polarizations of the neutron and proton den- 
sities of 208 Pb induced by Coulomb interaction. 



nuclei such as 208 Pb, this effect is weak while for the nuclei 
with N ~ Z this effect can be very distinguishable. 

We calculate the symmetry energy of 208 Pb where the 
Coulomb polarization effect of the nucleon densities is sub- 
tracted. The calculated results are listed in Table I. The ra- 
dial symmetry energy density distributions g(r) = 47^H sym r 2 
as a function of distance r are displayed in Fig. 2. Different 
interactions provide almost the same shape of the symmetry 
energy density distribution. The inset in Fig. 2 shows the nu- 
cleon density distribution in 208 Pb, from which one can see 
that the nuclear surface region is 6 - 8 fm away from the nu- 
clear centre. Though the matter in this region is much less 
dense than that in the nuclear interior, the contribution from 
the nuclear surface region to E sym (A) is rather large. The rea- 
son is twofold: Firstly, the spherical shell of the surface re- 
gion has a large volume. Secondly, the asymmetry distribu- 
tion function 5(r) = (p„(r) -p p (r))/p{r) has large values at the 
surface region. Fig. 3 shows the isospin asymmetry function 
6(r) versus distance r. Sir) is low in the interior of the nucleus 
but increases rapidly at the surface region. This also implies a 
neutron skin in 208 Pb. The region of r < 3 fm contributes a lit- 
tle to E sym (A) since it has only a small volume and S(r) is low 
in this region in spite of its high density. That the nuclear sur- 
face contributes greatly to a sym (A) indicate that the symmetry 
energy coefficient of a finite nucleus is considerably smaller 
than that of nuclear matter at saturation density due to the low 
density at the surface. 



TABLE I: The calculated symmetry energy coefficient a synl for Pb. 
The contributions of the gradient term a grad = C AifH smA r- dr and 
the second order term of the kinetic part a kin2 = f Q Anr 1( H^dr as 
well as their proportions in a sym are all listed, k is the ratio of the 
surface symmetry coefficient to the volume symmetry coefficient So- 





SIII 


SLy4 


SLy5 


SkM* 


a sym (MeV) 


21.6 


23.6 


23.3 


20.8 


fl grad (MeV) 


-0.290 


-0.348 


-0.399 


-0.307 


^grad f &sym 


-1.3% 


-1.5% 


-1.7% 


-1.5% 


flkin2 (MeV) 


0.365 


0.543 


0.516 


0.136 


#kin2 / Asym 


1.7% 


2.3% 


2.2% 


0.6% 


So (MeV) 


28.2 


32.0 


32.0 


30.0 


K 


1.81 


2.26 


2.21 


2.62 




polarizes the neutron and proton densities. One can calcu- 
late the neutron or proton density p by using the Skyrme in- 
teractions, and when the Coulomb interaction is switched off 
one obtains the density p . The differences between p and p , 
namely Ap„ = p„ —p n and Ap p — p p -p p , denote the Coulomb 
polarization effect. Fig. 1 displays the polarizations of the 
proton (upper panel) and neutron (lower panel) densities of 
2()8 Pb induced by the Coulomb interaction. The interactions 
used here (see in Fig.l) yield the binding energy for 208 Pb be- 
ing 1636.62 MeV, 1636.47 MeV, 1635.67 MeV and 1633.58 
MeV, respectively, which are in very good agreement with the 
experimental measurement of 1635.84 MeV 03611 . The protons 
and neutrons move slightly away from the core as a result of 
the Coulomb repulsion between charged protons. For heavy 



FIG. 2: (Color Online) Symmetry energy density distribution g(r) 



4^*H sym r 2 as a function of r in 208 Pb. 



The reduction of the a sym for 208 1 



Pb caused by the gradient 
term is less than 2%, as presented in Table I. Also, the second 
order term of the kinetic energy part contributes rather weakly 
to the a sym (less than 3%). Moreover, the contributions from 
these two weak terms tend to cancel out to a large extent. Thus 
one can neglect the two terms simultaneously without loss of 
accuracy. Accordingly, the Hj for finite nuclei may be sim- 
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to explore the density dependence of the nuclear symmetry 
energy coefficient S(p). The value of S(p) at p = 0.1 firT 3 
has been widely investigated recently. The analysis of the 
GDR of 208 Pb with Skyrme interactions suggests S(p = 0.10 
fm" 3 ) = 23.3 - 24.9 MeV JH while the analysis of the 
GDR of 132 Sn within relativistic mean field models leads to 
S{p = 0.10 fm" 3 ) = 21.2 - 22.5 MeV JH. Our calculated 
a sym (A) for 208 Pb is 20.8-23.6 MeV which agrees with that 
in Ref. (H] (20.22 - 24.74 MeV) obtained from the mea- 
sured nuclear masses, indicating the reliability of the present 
approach to some extent. 

III. DENSITY DEPENDENCE OF THE NUCLEAR 
MATTER SYMMETRY ENERGY FROM EXPERIMENTAL 
ALPHA-DECAY ENERGIES 



FIG. 3: (Color Online) Isospin asymmetry distribution function 
8(r) = (p„(r) - p p (r))lp(r) as a function of distance r in 208 Pb. 
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which leads to fisym = 'Ht + *Hv j ust a function of p(r) and 
S(r). For the large distance r, the asymmetry 6(r) is close to 
one, but this region almost does not contribute to the symme- 
try energy. Accordingly, this equation of the "Hj is valid. 

The mass dependence of the symmetry energy coefficients 
of nuclei is written as lUia 13711 
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where k is the ratio of the surface symmetry coefficient to 
the volume symmetry coefficient and Q is the surface stiff- 
ness that measures the resistance of a nucleus against separa- 
tion of neutrons from protons to form a neutron skin. Since 
both the a sym (A) and So can be estimated within the Skyrme 
interactions, the relation between the surface symmetry en- 
ergy and the volume symmetry energy characterized by the k 
value can be investigated directly. The k value of 208 Pb we 
obtained is 1.8 — 2.6, which agrees with that from the nu- 
clear masses [16] but slightly larger than that from the lep- 
todermous expansion based on the self-consistent mean-field 
theory [25] and from the binding energies in Ref. ll26ll . Re- 
cently, Centelles et al. proposed that the a sym (A) of finite nu- 
clei is approximately equal to S(pa) of the nuclear matter at 
a reference density pa, namely S(pa) = a sym (A) For 
example, the reference density for 208 Pb is p A = 0.1 fm~ 3 , 
namely, S(pa = O.lfm- 3 ) = a sym ( 208 Pb). This relation links 
the properties of finite nuclei and nuclear matter and helps us 



The nuclear symmetry energy coefficient a sym (A) is usually 
extracted by directly fitting the measured nuclear masses with 
different versions of the liquid-drop models Im |40[ EH . In 
this study, however, we extract a sym (A) with the experimental 
alpha decay energies of heavy (or superheavy) nuclei. In our 
previous work, a formula for Q a values was put forward which 
works well for heavy nuclei with Z > 92 and N > 140 I42I1 . 
That formula will be used here to extract a sym (A). And we 
use 162 heavy and superheavy nuclei in the extraction. The 
present approach is simple and physically clear since some 
less important terms in the mass formula of the liquid drop 
model are canceled out to a great extent in the calculation of 
Q a values. In this mass region, the shell correction energy 
is very small compared to the other terms [42]. So one is 
free from the influence of the shell effects. Subtracting the 
Coulomb term Qc, one obtains 

<N -Z\ 2 



Qa~Qc = ~- 



45 „ 



(id 



1 + kA- 1 ^ 

where Qc - aZA~ 4 t 3 (3A - Z), a — 0.9373 and e is a constant. 
The mass dependence of the symmetry energy coefficient is 
given by Eq. ( fTOt . 

To fit Sq and k directly is difficult to yield their optimal 
values since many different combinations of So and k could 
provide the same least deviation. Considering that S o has been 
relatively well determined as 31.6*2? MeV nowadays i2lTl . 
we solely determine the value of k in terms of a sym (A). The 
calculated value of k is 2.2 1 If w ^ a ^ ^ ne root-mean-square 
deviations of 0.33 MeV for Q a - Qc, where the uncertainty 
mainly results from the uncertainty of So value. The k value 
here from the experimental a-decay energies covers the range 
given by the Skyrme energy density functional method (see 
Table I). The two completely different approaches provide the 
unanimous results and hence they may validate each other to 
get more compelling results. 

Now let us constrain the density dependence of the nu- 
clear matter symmetry energy coefficient at subnormal den- 
sities with the help of the obtained a sym ( 208 Pb). We use the 
Eq. (6) in Ref. [21] to describe S(p), which reads 



/ \ 2/3 
S(p) = 17.47—1 +C 



Po) \P0 



1.52 



(12) 



5 



where C\ and C2 are linked by So = 17.47 + C\ + C2 and 
they need to be determined. We would point out that Eq. 
(fTZt is more universal for the description of the density de- 
pendence of Sip) 1I2TI1 compared with other forms. With 
S (pa) = o sym {A) and Eq. dT2b . we obtain 



1 +KA- 1 / 3 



17.47 



/ \ 2 / 3 



+ c 2 « 

\P0 



1.52 



(13) 



Inserting the detailed values of S 

fm- 3 for 208 Pb and k = 2.21+° ^, 
3+10.02 



31.6 +22 MeV,p A =0.10 



J -2.2 ' 

we obtain C\ = 24.56 



-7.82 
+7.58 



MeV and C 2 = -lOAltgji MeV, and accordingly the slope 
and curvature parameters are calculated to be L = 61 *™ MeV 



and /T sym 



-109!™ MeV. 



TABLE II: Values of the neutron skin thickness in 208 Pb obtained 
from various independent studies. 



Reference 



Method 



AR„ P (fm) 



Ref. f 151 pygmy dipole resonance 

Ref. f441 pygmy dipole resonance 

Ref. [451 proton inelastic scattering 

Ref. [46] proton elastic scattering 

Ref. [471 chiral effective field theory 

Ref. [48] exotic atom 

Ref. [49] electric dipole polarizability 

Ref. [211 mean-field 

Present alpha decay energies 



0.194 + 0.024 
0.180 ±0.035 

156 +0025 
211+ 054 

"■^ 11 -0.063 

0.17 ±0.03 
0.18 ±0.02 
0.168 ±0.022 
0.185 ±0.035 
0.191 ±0.032 



It has been well established in Ref. [43] that the correlation 
between L and AR np is AR„ P = 0.101 + 0.00147L for 208 Pb, 
and here L and AR„„ are measured in units of fm. With the 



values of L we obtained, the neutron skin thickness in 208 Pb is 
AR np = 0.191 + 0.032 fm in terms of this correlation. A com- 
parison of our result with that of others is presented in Table 
II. On the whole, our result is in agreement with the others, es- 
pecially the very recent data obtained from the pygmy dipole 
resonance. However, the values of the neutron skin thickness 
in 208 Pb obtained by various approaches carry large uncertain- 
ties currently. It seems that a measurement with higher accu- 
racy needs to be done in order to further well constrain the 
density dependence of the nuclear matter symmetry energy 
coefficient. 



IV. SUMMARY 

We have derived the density functional of the symmetry en- 
ergy for finite nuclei in the framework of the Skyrme-Hartree- 



Fock approach. The nuclear symmetry energy coefficient 
a sym (A) has been then directly extracted and investigated. The 
ratio of the surface symmetry coefficient to the volume sym- 
metry coefficient k was also estimated. For 208 Pb, a sym (A) was 
calculated to be 20.8-23.6 MeV with the Skyrme interactions 
SIII, SLy4, SLy5 and SkM*, and the corresponding ratio k is 
1.8 — 2.6. For a heavy nucleus, its surface region contributes 
largely to the symmetry energy according to the spatial dis- 
tribution of the symmetry energy. It has been shown that the 
contribution from the gradient term to the symmetry energy 
is found to be rather small (less than 2% for 208 Pb) and the 
second order term in the kinetic energy part of the symmetry 
energy density is also less important for the symmetry energy 
(its contribution is less than 3% for 208 Pb). Meanwhile, an 
alternative method (a macroscopic method) has been devel- 
oped in the present study to determine the symmetry energy 
coefficient of heavy nuclei which resorts to the available ex- 
perimental alpha decay energies of heavy nuclei. The value of 
k was calculated to be 2.21!°*? w j m $ = 31.6 + 2.2 MeV 
based on the macroscopic method, which agrees well with 
that from the Skyrme density functional approach. More- 
over, the two completely different approaches may validate 
each other to achieve more compelling results. The calcu- 
lated a sym (A) of 208 Pb with the macroscopic method was fur- 
thermore used to analyze the density dependence of the sym- 
metry energy coefficient of nuclear matter and evaluate the 
values of the slope and curvature parameters L and K sym . It 
was found that L = 61! 22 MeV and K sym = -109!™ MeV. 
The neutron skin thickness of 208 Pb was then estimated to be 
AR„ p — 0.191 + 0.032 fm, which is consistent with the pygmy 
dipole resonance studies. 



This work was supported by the Major State Basic Research 
Developing Program of China under No. 2013CB834405, the 
National Natural Science Foundation of China under Grants 
No. 11175219, 10975190, and 11275271; the Knowledge 
Innovation Project (KJCX2-EW-N01) of Chinese Academy 
of Sciences, CAS/SAFEA International Partnership Program 
for Creative Research Teams (CXTD-J2005-l)and the Funds 
for Creative Research Groups of China under Grant No. 
11021504. 



[1] P. Danielewicz, R. Lacey, and W. G. Lynch, Science 298, 1592 
(2002). 

[2] A. W. Steiner, M. Prakash, J. Lattimer, and P. J. Ellis, Phys. 

Rep. 411, 325 (2005). 
[3] V. Baran, M. Colonna, V. Greco, and M. Di Toro, Phys. Rep. 



410, 335 (2005). 

[4] B. A. Li, L. W. Chen, and C. M. Ko, Phys. Rep. 464, 1 13 (2008). 
[5] J. M. Lattimer and M. Prakash, Phys. Rep. 442, 109 (2007). 
[6] Jianmin Dong, Wei Zuo, and Werner Scheid, Phys. Rev. Lett. 
107, 012501 (2011). 



6 



[7] C. J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 86, 5647 
(2001). 

[8] Bharat K. Sharma, Subrata Pal, Phys. Lett. B682, 23 (2009). 
[9] J. M. Lattimer and M. Prakash, Phys. Rep. 333, 121 (2000); 
Science 304, 536 (2004). 
[10] B. G. Todd-Rutel and J. Piekarewicz, Phys. Rev. Lett. 95, 
122501 (2005). 

[11] C. J. Horowitz, S. J. Pollock, P. A. Souder, and R. Michaels, 

Phys. Rev. C 63, 025501 (2001). 
[12] T. Sil, M. Centelles, X. Vinas, and J. Piekarewicz, Phys. Rev. C 

71, 045502 (2005). 
[13] Isaac Vidana and Constanca Providencia, Artur Polls, Arnau 

Rios, Phys. Rev. C 80, 045806 (2009). 
[14] Lie-Wen Chen, Che Ming Ko, and Bao-An Li, Phys. Rev. Lett. 

94, 032701 (2005); Phys. Rev. C 72, 064309 (2005); Bao-An 

Li and Lie-Wen Chen, Phys. Rev. C 72, 064611 (2005). 
[15] Andrea Carbone et al., Phys. Rev. C 81, 041301(R) (2010). 
[16] M. Liu, N. Wang, Z. X. Li, and F. S. Zhang, Phys. Rev. C 82, 

064306 (2010). 

[17] M. B. Tsang et al., Phys. Rev. Lett. 102, 122701 (2009). 
[18] P. Danielewicz and J. Lee, Nucl. Phys. A818, 36 (2009). 
[19] M. Centelles, X. Roca-Maza, X. Vinas, and M. Warda, Phys. 

Rev. Lett. 102, 122502 (2009). 
[20] M. Warda, X. Vinas,X. Roca-Maza, and M. Centelles, Phys. 

Rev. C 80, 024316(2009). 
[21] Jianmin Dong, Wei Zuo, Jianzhong Gu, and Umberto Lom- 

bardo, Phys. Rev. C 85 034308 (2012). 
[22] B. A. Brown, Phys. Rev. Lett. 85, 5296 (2000). 
[23] S. Typel and B. A. Brown, Phys. Rev. C 64, 027302 (2001). 
[24] R. J. Furnstahl, Nucl. Phys. A706, 85 (2002). 
[25] P.-G. Reinhard, M. Bender, W. Nazarewicz, and T. Vertse, 

Phys. Rev. C 73, 014309 (2006); N. Nikolov, N. Schunck, W. 

Nazarewicz, M. Bender, and J. Pei, Phys. Rev. C 83, 034305 

(2011). 

[26] A. E. L. Dieperink, and P. Van Isacker, Eur. Phys. J. A 32, 11 
(2007). 

[27] Kazuhiro Oyamatsu and Kei Iida, Phys. Rev. C 81, 054302 
(2010); Kazuhiro Oyamatsu, Kei Iida, and Hiroyuki Koura, 
Phys. Rev. C 82, 027301 (2010). 

[28] Ning Wang and Min Liu, Phys. Rev. C 81, 067302 (2010). 



[29] K. Bennaceur, J. Dobaczewski, Comput. Phys. Commun, 168, 
96 (2005); M.V Stoitsov, J. Dobaczewski, W. Nazarewicz, P. 
Ring, Comput. Phys. Commun, 167, 43 (2005), and references 
therein. 

[30] S. Cwiok, P.-H. Heenen, and W. Nazarewicz, Nature (London) 

433, 705 (2005), and references therein. 
[31] J. R. Stone, P.-G. Reinhard, Prog. Part. Nucl. Phys. 58, 587 

(2007) , and references therein. 

[32] S. Goriely, N. Chamel, and J. M. Pearson, Phys. Rev. Lett. 102, 
152503 (2009); Phys. Rev. C 82, 035804 (2010), and references 
therein. 

[33] J. Bartel and K. Bencheikh, Eur. Phys. J. A 14, 179 (2002). 
[34] M. Brack, B. K. Jennings and Y.H. Chu, Phys. Lett. B65, 1 
(1976). 

[35] B. Grammaticos and A. Voros, Ann. Phys. 129, 153 (1980). 
[36] P. Moller, J. R. Nix and K.-L. Kratz, At. Data Nucl. Data Tables, 
66, 131 (1997). 

[37] W. D. Myers and W. J. Swiatecki, Ann. Phys. (N.Y.) 55, 395 

(1969); 84, 186 (1974). 
[38] L. Trippa, G. Colo, and E. Vigezzi, Phys. Rev. C 77, 061304(R) 

(2008) . 

[39] L. G. Cao and Z. Y. Ma, Chin. Phys. Lett. 25, 1625 (2008). 
[40] Ning Wang and Min Liu, Phys. Rev. C 81, 067302 (2010). 
[41] T. Mukhopadhyay and D. N. Basu, Acta Phys. Pol. B38, 3225 
(2007). 

[42] Jianmin Dong, Wei Zuo, Jianzhong Gu, Yanzhao Wang, and 

Bangbao Peng, Phys. Rev. C 81, 064309 (2010). 
[43] X. Roca-Maza, M. Centelles, X. Vinas, and M. Warda, Phys. 

Rev. Lett. 106, 252501 (2011). 
[44] A. Klimkiewicz et al., Phys. Rev. C 76, 051603(R) (2007). 
[45] A. Tamii et al., Phys. Rev. Lett. 107, 062502 (2011). 
[46] J. Zenihiro et al., Phys. Rev. C 82, 04461 1 (2010). 
[47] K. Hebeler, J. M. Lattimer, C. J. Pethick, and A. Schwenk, 

Phys. Rev. Lett. 105, 161102 (2010). 
[48] E. Friedman and A. Gal, Phys. Rep. 452, 89 (2007). 
[49] J. Piekarewicz, B. K. Agrawal, G. Colo, W. Nazarewicz, N. 

Paar, P.-G. Reinhard, X. Roca-Maza, and D. Vretenar, Phys. 

Rev. C 85, 041302(R) (2012). 



